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ABSTRACT 



Initial design study of a low altitude satellite (LASAT) communication sys- 
tem is performed. The use of many inexpensive, low-altitude satellites in random 
orbits may be a solution to vulnerability of current military communication sys- 
tems. Statistical study of orbit characteristics is performed giving the mean 
number of satellites in view of a ground station and the coverage density of the 
satellite array. Waveform analysis is performed on a coded and uncoded, 
orthogonal, noncoherent, fast-frequency-hopped M-ary frequency-shift-keyed 
signal in a Rician fading channel with optimum partial band jamming. An ana- 
lytical expression for the system probability of bit error is obtained and numerical 
results are generated for various lex els of fading, jamming, and diversity. For- 
ward error correction coding is applied using convolutional codes and Reed- 
Solomon codes. 
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I. INTRODUCTION 



Current U. S. military communciations are performed using a small number 
of highly sophisticated and expensive satellites in geostationary orbits. The loss 
of only a few of these satellites will have a serious impact on vital communi- 
cations. Although current technology has not yet produced an anti-satellite 
(ASAT) system that can reach geostationary orbits, such systems may be devel- 
oped in the future. In addition, ground tracking and control stations are suscep- 
tible to breakdown and attack. To close this vulnerability window many more 
satellites must be fielded quickly. The problem is the enormous cost of building 
and launching these high-altitude satellites, which is estimated at hundreds of 
millions of dollars each. 

A possible solution to this vulnerability in times of crisis is the use of less so- 
phisticated. low-altitude satellites designed to be launched quickly in large num- 
bers. These satellites would be light and inexpensive to build and operate. 
Hundreds of these satellites would be placed in orbit to provide full earth cover- 
age and their sheer numbers would overwhelm any ASAT system. To achieve 
reasonable costs, the satellites will have minimum station-keeping ability, if any. 
Initial system analysis has indicated that about 250 satellites in random orbits 
would be needed to provide full earth coverage with inplace reconstitution 
[Ref. 1]. The multi-satellite system (MSS) considered in Ref. 1 is primarly ex- 
pected to be used as a strategic replacement to the current geostationary system 
and is designed to pro\ide world-wide communications. This configuration 
would require that the satellites communicate between each other via a packet 
switching system. The system considered in this thesis is similar to the multi- 
satellite system of Ref. 1 but would be used for short term tactical communi- 
cations. Therefore, the satellite crosslinks are unnecessarv. 



II. BACKGROUND INFORMATION 



A. LOW-ALTITUDE SATELLITES 

The main purpose of Low-Altitude Satellites (LASAT's) is to handle crisis 
communications for military commanders. The satellite is. expected to operate 
with encrypted data and digitized voice links to local areas, possibly between 
hand-held computers in the field. One advantage of LASAT is the small propa- 
gation delay as compared to 250 ms for geostationary satellites. This will facili- 
tate the speedy transmission of data, and many data protocols designed for 
terrestial radio communications can be modified for LASAT use. 

On the other hand. LASAT coverage is very small, and the window of 
communciations will last for only minutes instead of hours. This requires a great 
number of satellites to provide continuous coverage. It is envisioned that 
LASAT's will be deployed randomly in many orbits without station-keeping ca- 
pability and that communications will be through omnidirectional antennas. 
Such a random system requires more satellites to provide continous coverage than 
does one with well-defined orbits where station keeping is required. Therefore, 
the number of satellites seen by a ground terminal is a random variable, and a 
ground station signal may be received simultaneously by many satellites. Also, 
when used for tactical communications without cross-links, the same ground 
signals relayed by many LASAT's become satellite multipath signals at the re- 
ceiving ground terminals with near identical power if classical transponders are 
used. In addition, because low-gain antennas will be used by the ground termi- 
nals and LASAT's, reflection multipath signals are also present. The combina- 
tion of these multipath effects result in a fading channel. 

B. THE FADING CHANNEL 

Generally, radio waves do not take a direct path from the transmitter to the 
receiver. Terrain, buildings, vegetation, and the ionosphere cause multiple re- 
fections and refractions of the signal. Therefore, a communications channel is. 



in most cases, a multipath channel with the signals arriving at the receiver from 
a variety of directions and with different delays. Furtheremore, the channel is 
very dynamic with the movement of the ionosphere, receiver and transmitters, 
and foliage causing a time variance in the structure of the medium. As a result 
of such time variations, the multipath varies with time. Small changes in path 
length can also result in large phase changes in the signal at the receiver antenna. 

Stein [Ref. 2] describes the system in phasor terms, where the observed re- 
ceived phasor is a vector sum of several phasors with each varing randomly and 
individually over a full 2 r, range. As a result, constructive and destructive inter- 
ference occurs and the received signals will have different amplitudes. The fading 
channel, therefore, is characterized by a random time-variant impulse response. 
Thus, the transmission of a pulse will result in a pulse train arriving at the re- 
ceiver with variations in amplitude. If a series of pulses are sent, the pulse train 
as* elated with each transmission will be different as shown in Figure 1 on page 
4 [Ref. 3 : p. 704]. These time variations are unpredictable to the user and may 
be considered random which implies that the received signal can be modeled as 
a random process. 

In the general case, the received signal will have both a direct ( or specular) 
component and a number of scattered (or diffuse) components. The scattered 
components' arrival times and amplitudes will vary randomly about the mean 
provided by the direct component. As a result, the envelope of the received signal 
will be similar to the envelope of a non-zero mean signal. This type of channel 
is described as Rician. The amplitude of the received signal has the probability 
density function (pdf) [Ref. 4: p. 105] 




where cr represents the power of the direct component, 2 o'- is the expected value 
of the power of the scattered component, and I 0 (x) is the zeroth-order modified 
Bessel function of the first kind. For the case where there is no direct component 
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figure 1. Example of a multipath channel on transmitted pulses: from 

Ref. 3: p.704. 



and all the received energy is via the scattered components, the channel is de- 
scribed as Rayleigh with the pdf of the amplitude given by 



/w(< ? ) = 




a > 0. 



(2.2) 



The fading characteristic of the channel may be rapid or slow dependent 
upon the relative time and frequency scale of the medium and the signal band- 
width. The coherence bandwidth is defined as a measure of the frequency 
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dependence of the channel. If the channel coherence bandwidth is smaller than 
the bandwidth of the transmitted signal, the channel is frequency-selective and 
the signal is severely distorted. However, if the coherence bandwidth is larger 
than the signaling bandwidth, the channel is frequency-nonselective, and the sig- 
nal is undistorted. A slowly-fading channel is one in which the coherence time is 
large and the signal is unchanged over a signaling interval. On the other hand, 
in a rapidly-fading channel, the signal is distorted by changes in amplitude or 
frequency during the signaling interval. 

Unusual conditions in the ionosphere can result in significant scintillations 
of earth satellite links. These effects are highly frequency dependent. In mild 
scintillation the channel is typically Rician but at times of strong scintillation 
(deep fading) has been observed to be Rayleigh. The amount of scintillation is 
dependent upon satellite altitude, time of day, and latitude, with the largest ef- 
fects near the nighttime equatorial zone. Further significant multipath effects 
that would be involved in a LASAT system are dependent upon terrain and the 
speed of the satellite and terminals. The relative speed between the satellite and 
the ground station can be significant with a low earth orbit satellite system, and 
will increase fading effects. [Ref. 2] 

Diversity is one technique that can be used to overcome the effects of fading. 
With diversity the same information is provided over two or more independently 
fading channels. This lowers the probability that any individual symbol will be 
affected by deep fading on each channel. There are a variety of methods of di- 
versity combining and they include the following [Ref. 2]: 

® Spaced antennas 

Antenna spacing of ten wavelengths or greater can result in independent 
channels. This technique, however, is not well suited to satellite UHF com- 
munications because of the long wavelengths. 

® Frequency diversity 

The use of wideband frequency hopping (fast hopping) on fading chan- 
nels. with the same information carried on different hops, tends to result in 
independent states of fading. 



• Multipath diversity 

A signal having a bandwidth greater than the coherence bandwidth of 
the channel will resolve multipath components and thus, provide the the re- 
ceiver with several independent fading signal paths. This is an inherent 
characteristic of direct sequence spread spectrum systems. 

• Time diversity 

The same data is repetitively transmitted at intervals that exceed the 
coherence time of the channel which results in independent fading condi- 
tions. 

Coding and interleaving are other techniques which are similar to time di- 
versity in providing redundancy. This is particularly true with binary linear 
block codes. However, convolutional codes and other non-binarv linear codes 
provide satisfactory results. A more detailed coverage of fading channels is pre- 
sented in Ref. 2 and 3. 

C. FREQUENCY HOPPING SPREAD SPECTRUM 

The previous section discussed the effects of the fading channel and indicated 
some methods available to overcome the resulting degradation in system per- 
formance. Fast frequency hopping spread spectrum, with the same information 
carried on different hops, is one form of frequency diversity that may be used. 
If sequential hops exceed the coherence bandwidth of the channel, then fading for 
each hop frequency is independent. The selection of the frequency slots in each 
signaling interval is made pseudo-random according to the output of a pseudo- 
noise generator. 

Due to rapid signal phase changes which occur within a fading channel and 
the difficulty in building and using coherent demodulators for frequency shift 
keyed systems, most fast frequency hopping systems use noncoherent, orthogonal 
modulation schemes. A common data modulation for FH systems is M-ary fre- 
quency shift keying (MFSK). were K = log ; M bits are used to determine which 
of M frequency tones will be used for carrier modulation. The tones are usually 
spaced far enough apart so that the transmitted signals arc orthogonal. This 
implies that the frequency spacing be at least 1 j KT b Hz. where T b is the bit du- 
ration and T, = KT b is the duration of each symbol. 
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When fast frequency hopping is used with a M-ary FSK system, each of the 
M symbols is subdivided into L chips. At each chip, the MFSK modulator out- 
put is hopped to a different frequency. Since the chip duration T c is shorter than 
the data modulator output symbol duration T r the minimum tone spacing for 
orthogonal signals is B = 1 / T c = L/ KT b . The total system bandwidth then be- 
comes IV = m.MB = mLM/KT b , where m is an arbitary integer that defines the 
hopping bandwidth. 

A typical frequency hopping system is shown in Figure 2 on page 8. At the 
receiver, an identical PN generator synchronized with the received signal removes 
the pseudo-random frequency translation by mixing the synthesizer output with 
the received signal [Ref. 5]. For an orthogonal, noncoherent system, demodu- 
lation is usually accomplished through square-law, or envelope detectors. The 
detector output can then make a hard decision at each chip interval as to which 
of the M-channels contains the signal and make an estimate on all L chip deci- 
sions, or it can sum the energy of L chips on each channel and make a soft deci- 
sion based on the largest value. 

Although the fast frequency hopping system is effective in a fading channel, 
partial band jamming can seriously degrade the signal when it hops in and out 
of the jamming band. The partial band jammer concentrates its power over a 
small fraction ?/ of the system bandwidth. On some hops, the signal will be at a 
frequency which is jammed, while on other hops, the signal will be at a frequency 
where the only interference is thermal noise. A system which can measure the 
amount of thermal and jamming noise power per hop can provide this side in- 
formation to normalize the measurements obtained at the detector outputs 
(chapter IV). Worst case jamming occurs when the jammer optimizes 7/ to maxi- 
mize the system bit error probability. In system design it is normally assumed 
that the jammer has knowledge of the system state and is able to optimize ?/ . 
Therefore, performance calculations provide an upper bound on error rates. 

D. SUMMARY OF RESEARCH 

This thesis covers the communications analysis of LASAT system design. 
Included in this studv arc orbit characteristics of the svstem. and the waveform 
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Figure 2. Typical Frequency hopped spread spectrum system 

analysis of a proposed fast frequency hopped M-ary frequency shift keyed mod- 
ulation scheme. 

1. Orbit statistics 

Chapter 111 provides a study of the LAS AT orbit distribution and cover- 
age. The first section of the chapter looks into the number of satellites in view 
of a ground station. An analytical approach was used to derive the mean value 
and the standard deviation of the number in view. A computer simulation veri- 
fied the results. The second section of chapter III discusses the coverage of the 
LASAT relative to a ground terminal. The mean coverage area and the standard 



S 





deviation in which connectivity can be established with other ground terminals 
was found. A computer simulation verified the results. 

2. Waveform selection 

Waveform selection involves spread spectrum techniques, carrier modu- 
lation. error correction coding, and diversity. Chapter IV presents the derivation 
of the bit error probability for uncoded and coded fast-frequency-hopping M-ary 
frequenev-shift-keying schemes in a Rician fading channel with optimum partial 
band jamming. Numerical results for various levels of diversity and fading arc 
provided for both the coded and uncoded systems. Convolutional codes and 
Reed-Solomon codes were used to provide forward error correction coding. 
Chapter V presents conclusions and recommendations for further research. 
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III. LASAT SYSTEM 



A. DISTRIBUTION OF LASAT 

Since the LASAT communication system is envisioned as consisting of a large 
number of low earth orbiting satellites distributed randomly throughout many 
orbits, it is expected that the number of satellites in range of a ground station will 
also be random. The distribution of satellites within view of a ground station will 
be developed here both theoretically and numerically. The ground station is as- 
sumed to have an omnidirectional antenna that will reach satellites with an ele- 
vation angle of zero degrees. A general structure of L satellites randomly 
distributed within an orbit altitude range /?, to h 2 will be assumed, 
i. Theoretical development 

The total volume J ' enclosed by the orbital radii R , and R 2 is 

r = Y^2-^) (3.i) 



where 



R i = R e + h \ ,, 

R 2 = R e + h 2 

with R e = 6378.14 km. (3443.93 NM). the radius of the earth. The volume den- 
sity of L satellites within this volume is 

p — Ljy satellites 'km 3 . (3.3) 

The volume of the orbit portion seen by a ground station can be derived with the 
help of Figure 3 on page 12. The volume of the spherical sector v sec2 defined by 
points ACBO with radius R : is 
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where 



cos u = R e j{R e + h 2 )> 



(3.5) 



The right circular cone defined by the points ABO has a volume of 



' cone 



y V 

y - Rc)K- 



(3.6) 



where a is the radius of the base and a 2 = Rt — R; . The volume of the segment 
and base ABC is the difference between the volumes of the spherical sector and 
the right circular cone. 



' 2 ' scc2 ' cone 

= y ~{R e + h 2 f - y ~[(y + luf - R^R e 
= y r[2 R‘h 2 + ARM; + 2 h] - ( R ~ + 2 RJu + h 2 R e - /$] ^ 

— - 7 - n/z-, [3/?-. — 3/C A 2/A] 

J) " - 

= 4 - r.h;(2>R, - Ik). 

By a similar method, the portion of volume seen by the ground station below /?, 
is 
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Figure 3. Geometry of satellite orbit distribution. 



= ~r~ /'| 2 (3/?, — /i,), (3.8) 

and the subvolume v of the orbital shell of thickness h 2 —h\ seen by the ground 
station is 
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(3.9) 



Therefore, on average, the number of satellites seen by the ground station will be 
the density of satellites times the subvolume v, or by using (3.3) and (3.9) 



The above result can also be shown statistically with the help of the de- 
velopment by Larson [Ref. 6: pp. 143-149]. The distribution of satellites 
throughout the volume V is random and uniform and one can assume that a 
volume v can be divided into m non-overlapping subvolumes of size Av. where 
Av = vjm is very small. It is therefore reasonable to assume that there will or will 
not be a satellite in each of these small volumes and that the probability that ex- 
actly one satellite will be in each subvolume is pAv = pvjm for each, where p is 
the density of L satellites distributed throughout the total volume V. Further- 
more. it can be assumed that the presence of a satellite in each subvolume is in- 
dependent of the other subvolumes. Thus the volume v has been subdivided into 
m repeated independent Bernoulli trials, each of a volume Av = vjm with the 
probability of a satellite being present equal to pAv = pv/m. If n is the number 
of satellites in the subvolume v. then n is approximately a binomial random vari- 
able with probabiltity density function /%(«) where 



n 



ave 




(3.10) 




n = 0.1.2 m 



(3.11) 



= 0. otherwise. 
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If the expression is taken to the limit m — * oo (thus the individual subvolumes are 
shrinking to zero), the exact probability function for the Poisson random variable 
results. That is 



*<«>-&: x-£)o~£) 



m—n 



{py) n 

n\ 



-pv 



(3.12) 



Therefore, the distribution of the number of satellites within a volume v can be 
represented by the volumetric Poisson distribution with a parameter pv. The 
well-known mean n and variance c : n of the Poisson process is 

n = pv = -y- , (3.13) 



nd 







= pv = 




(3.14) 



2. Numerical simulation 

To generate a sample set representative of the distribution over the orbital 
volume, a uniform random sequence of L points between 0 and I’ is generated. 
To find the number of satellites seen by the ground station for a particular sample 
set of L points between 0 and J\ the number of points between 0 and v are 
counted, where V and v the orbital volume and subvolume defined by (3.1) and 
(3.9). When a large number of sample sets are run and the results averaged, the 
mean number of satellites n and the variance o'- of the samples can be found from 

M 

( 3 - 15 ) 

;=1 



and 
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where M is the number of sample sets and n,- is the number of satellites found 
within the subvolume v in the i th sample set. 

Sample means were found by generating first 2000, and then 5000 sample 
sets of 250 and 150 points uniformily distributed as described above. The num- 
ber of satellites found in the sub-area were then averaged and compared against 
the mean found by equation (3.13). Table 1 shows the numerical results for 
/z, = 550 km (300 NM) and lu = 750 km (400 NM) with 2000 and 5000 samples. 
The numerical results compare favorably with the theoretical results obtained by 
(3.15) and (3.16). Differences are due to the pseudo-randomness of the uniform 
number generator used. 



Table 1. MEAN AND VARIANCE OF SATELLITE DISTRIBUTION. 





Theory 


2000 Samples 


5000 Samples 


L 


n 


cl 


n 


cl 


n 


* s 

c;, 


250 


1 1 .54 


1 1 .54 


12.51 


10.43 


12.4S 


10.89 


150 


6.92 


6.92 


7.89 


6.78 


7.87 


6.52 



When an experiment described by a Poisson process is repeated a great 
number M of times, then it can be expected that the number M n of times that 
exactly n satellites are observed will be [Ref. 7: p. 158] 

U r *Mf s {n). (3.17) 

To verify this, the number of satellites found within the subvolume for the run 
of M = 5000 sample sets was plotted versus the number of times (occurences) M n 
that a particular number of satellites were within view (see Figure 4 on page 



17). This curve is compared to a Poisson distribution with the parameter pv as 
given by (3.13), and shown in Table 1 (see Figure 4 on page 17). 

Table 2 shows the maximum number of satellites seen by the ground 
station with a .999 and .990 probability. Since the distribution is Poisson, the 
probability P that there will be no greater than satellites seen by the ground 
station is 



P = Pr[/z < /t max l 

/2 max 

VW" - 



pv 



n\ 



77=0 



(3.18) 



The probability P 0 that a ground station sees no satellites overhead is (3.18) with 
'W = 0. or 

Po = e- p \ (3.19) 



and is also given in Table 2 



Table 2. NUMBER OF SATELLITES SEEN BY GROUND STATION 



L 


Hnxx 


^max 


n = 0 


P = 0.999 


P = 0.99 


Po 


250 


22 


20 


0.00001 


150 


16 


14 


0.001 



The above analysis shows the expected number of satellites a ground sta- 
tion may communicate with. For a transmitting receiving pair of ground 
stations, the number of satellites that may be used for communicating is depend- 
ent upon a joint pdf between the number of satellites, the orbital altitude and the 
distance between the stations. An analysis of the satellite coverage will be cov- 
ered in the next section. 
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Figure 4. Distribution of the number of satellites within view: Compared with 

ideal Poisson distribution with parameter pv= 11.54 . 



B. COVERAGE 



1. Theoretical development 

Now that the mean number of satellites seen by a ground station has been 
established, it is desirable to develope a distribution for the coverage area of the 
satellites. A satellite on the horizon of the base station has a sub-point that is an 
angle a away from the base station (sec Figure 5 on page 19), where 



The satellite at point S can communicate an additional distance of arc length 
R c y. to a user terminal at point UT. Therefore, the base station at point BS has 
a maximum communication range of arc length 2 c/.R e out to the user terminal 
through the satellite. This gives a maximum coverage area defined by the solid 
angle of 4y. or 4 cos ~ ] (RJR 2 ) steradians, which for an average altitude of 
h = (h : + /?[)/ 2 = 650 km (350 NM) is an area of 



Figure 6 on page 20 shows an example of the coverage area provided by 
a number of satellites. The spherical caps of area covered by the satellites are 
shown as circles of an equivclent area. In this diagram there are three satellites. 
SI. S2. and S3, uniformly spaced over the total view area of the base station 
(spherical cap of a radius of arc length a.R e from the base station). Each satellite 
has a range of arc length v.R e also. User terminals. Tl, T2. T3, T4. and T5, are 
also randomly spaced ox er an area out to a range of arc length 2 y.R c from the base 
station. One can sec from this Figure 6 that the user terminals hax'e differing 
numbers of satellites through which they can communicate to the base station. 
In this example, the user terminals can communicate as folloxvs: T4 through three 
satellites. T2 through two satellites. Tl and T5 through one satellite each, and 
T3 has no satellite in range. In general, as the range between the user terminal 
and the base station increases, there are fewer satellites through which to 



a = cos \RJR 2 ). 



(3.20) 



A = R; cos~' iRJ[R c + !,)■] 
= 1 7.63 x 10 6 km". 



(3.21) 
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Figure 5. Geometry of satellite coverage. 

communicate. If more satellites were added to the system, the communication 
ablilitv would increase. 

The coverage density for the total number of satellites L in the orbital 
volume V needs to be determined. This density represents the number of satel- 
lites that can communicate with both a user terminal and the base station as a 
function of the distance r between them for a given number n of satellites that arc 
within view of the base station. The satellites arc uniformly distributed through 
the orbital volume and their different altitudes will result in coverage areas of 
different sizes. This adds another random variable to the process. Since the 



19 




altitude difference is assumed to be small compared with the coverage area of 
each satellite, to simplify calculations without greatly affecting the results, it is 
safe to assume that all the satellites are at the same altitude and that the orbits 
are circular. The mean altitude h will be used in calculations. The area covered 
by the satellites can be modeled as an X-Y cartesian grid with the base station 
at the origin. The positions of the individual satellites arc represented by coor- 
dinate points Xj and y„ which are uniformity distributed from — c/.R c to a R e along 
the X and Y axes. Therefore, the pdfs for the random variables x and y are 

f x {x) = 1 !2uR e , - a.R e < x < aR e 

fy{y ) = 1 l2v.R e , - uR e < y < a R e . U j 

with a zero mean and a second moment E[.v] = E[y] = o 2 x = aj. - (aR e y/2. 

For a given number of satellites n, there are n sets of randomly distributed 
coordinates describing the satellite positions, x u y,.i= 1,2 . Each of the 
satellites covers an area defined by the radius of arc length c/.R e around the 
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satellite subpoint. The average distance from the base station to a satellite 
subpoint can be described by two new random variables 



n 




IY* 



n 




where z x and z y are the average X and Y axes coordinates for the satellites. Since 
n is a Poisson random variable, the distribution of z x and z y will be conditioned 
on /?. The conditional probability density functions of z x and z y are given by 



where the * is the convolution function. The central limit theory provides that if 
n is sufficiently large (greater than five), the sum of the random variables ap- 
proach a Gaussian distribution. Thus. z x and z y can be approximated as inde- 
pendent indenticaliy distributed Gaussian random variables with zero mean and 
a variance of c] 




( 3 . 24 ) 




( 3 . 25 ) 




where 
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(3.26) 



o] = E[zl] = E [zh 



= E 



XJ L }' 

n n 



-r ES a >'> = -T ZZ E M 

n /=1 j=l J n i—l j—\ 



£I.v-j = 



n 



w 

3 n 



i=i 



Since z T and z are independent, their joint density function is given by 



fz x z y I ,Mx' Z y I n ) = fz x |«( Z U n Y Zy I nYy I «) 



2rr(j 



i f z ; + 



(3.27) 



The base station can communicate with an area in the shape of a circle 
with radius 2 y.R e and the distribution can be represented as a function of the 
distance r of the user terminal to the base station. To do this, define a new ran- 
dom variable 



(3.2S) 



which is the range from the base station to the user terminal, and 



F R\Y'\n) - f 7 ^ | n {z x ,z y I n)dz x dz y 



[: x : y )sC 



2~o~. J 



exp' - 



+ z y 






2(7 2 . f 






(3.29) 



The region C is a circle of radius r. The above integral can be easily solved using 
polar coordinates [Ref. S: p. 9j. Let 

z x = r cos 0 
z v = r sin 0 
dz x dz y — rdrdO. 



*>0 



(3.30) 



